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ABSTRACT 


A near miss underwater explosion can significantly damage 
improperly shock hardened combat systems equipment and render 
the ship uncible to “fight hurt*'. MIL-S-901D currently requires 
shock qualifying mediximweight equipment to a "generic" shock 
excitation on the Navy's Mediumweight Shock Machine (MWSM). 
This shock excitation is severe, but not always characteristic 
of the actual ship structure response to an underwater 
explosion. This study proposes a design modification which 
will allow using a multi-DOF equipment mounting fixture on the 
MWSM which can be "tuned" to simulate shipboard shock 
characteristics determined from modal testing or previous ship 
shock trial data. Equipment qualified in this manner could be 
highly relied on to survive in battle. 
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1 . 


INTRODUCTION 


In today's world of modern warfare the U.S. Navy must rely 
on state of the art combat systems equipment to maintain its 
fighting advantage over possible advisaries. However, the 
high-tech nature of this equipment can lead to increased 
vulnerability to mechanical shock induced failure if not 
adequately packaged to withstand the severe excitations 
expected from conventional and nuclear underwater weapons. 
When exploded in proximity of the ship, these weapons produce 
an intense pressure wave v;hich is applied over the entire 
underwater portion of the ship's hull. Although the hull is 
likely to remain intact, the violent complex shock waves that 
propagate throughout the ship can sufficiently damage 
essential equipment to render the ship unable to "fight hurt" 
and significantly inpair its mission integrity. 

The Navy currently uses Military Specifications 
(MIL-S-901D), “Shock Tests, High Impact; Shipboard Machinery, 
Equipment 6uid Systems, Requirements For" to detail specific 
shock qualification reejuirements for shipboard machinery, 
equipment, systems and structures which are required to resist 
the effects of mechanical shock. These requireme.nts establish 
the general shock test criteria and provide the contracting 
activity a basis for selecting the appropriate testing device 
based on the weight category of the equipment. The weight 
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rjategories defined are; lightweight for an attached weight up 
to 550 lb; mediumweight for an attached weight up to 7400 lb; 
and heavyweight for a total weight up to 60,000 lb. The 
lightweight and mediximweight machines are similar in that the 
high impact, shock is delivered to the attached equipment by 
use of a hammer and amvil assembly. Hea'/yweight category 
equipment is installed onboard a floating platform barge and 
subjected to shock from an underwater explosive of known 
charge and standoff geometry. Of the three weight categories 
listed, this research focuses on the use of the Medium Weight 
Shock Machine (MWSM) shown in Figure 1. 



Figure 1. The Navy High-Iitpact Shock Machine for 
Mediumweight Equipment (MWSM). Courtesy of Clements(1972). 
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The MWSM delivers a vertical high impact mechanical shock 
to the anvil table by means of a 3000 lb swinging hammer. The 
impact between the hammer auid anvil is highly elastic and the 
energy is controlled by adjusting the height of the hammer 
above the anvil table prior to release. Required hammer 
heights are specified in MIL-S-901 based on the total weight 
attached to the anvil table. 

Deck mounted surface ship equipment is normally attached 
to the anvil table using the standard mounting fixture which 
is designed to provide a great deal of flexibility in 
equipment mounting geometry. The standard fixture is shown in 
Figure 2. Submarine deck mounted equipment is attached using 
the coil spring soft deck simulator shown in Figure 3. When 
using the standard fixture, MIL-S-901D specifies the number 
and type of support channels to be used based on the equipment 
weight and mounting bolt spacing. As noted by Clements(1972), 
the arrangement specified in MIL-S-901D was design to keep the 
calculated meiximum stress in the cheuinels below 35,000 psi in 
a static acceleration field of 50 g*s. Although not by design, 
using this arrangement produces an equipment excitation in the 
range of 55 Hz to 72 Hz. The soft deck simulator has a natural 
frequency in the reinge of 20 Hz to 25 Hz and is designed to 
simulate the natural deck frequencies of a submarine. 

Equipment tested using the standard fixture will be 
excited predominately at a single frequency in the 55 Hz to 72 
Hz range which is not realistic of the excitation aboard ship. 
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Figure 3. Soft Deck Simulator for MWSM. Courtesy of Hughes 
Aircraft. 
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Onboard ship, the excitation will be a coit^lex combination of 
the natural frequencies of the ship's structure in the 
location of the equipment installation. Shipboard vibration 
may excite natural frequencies within the equipment and lead 
to severe damage that would have not been detected on the 
MWSM. This introduces the need to develop a replacement for 
the standard fixture which will produce a shock excitation 
more realistic of the local shipboard response to underwater 
explosions. 

This study examines the design of a multi-degree-of- 
freedom fixture to be installed on the MWSM which can be tuned 
to more accurately simulate shipboard shock excitation of 
combat systems equipment. The shock spectrum of the excitation 
input to three representative pieces of equipment modeled in 
the DDG-51 Class Ship Pre-Shock Trial Analysis performed by 
Costanzo and Murray (1991) at the Underwater Research Division 
of David Taylor Research Center (DTRC/UERD) is used to 
identify the desired excitation response of the fixture. These 
shock spectra identify the frecjuency components of the 
equipment excitation eind provide the basis for tuning the 
fixture. This fixture is proposed as an economical means of 
providing realistic excitations to combat systems equipments 
when shock qpialifying on the U.S. Navy's MWSM. 
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II. BACKGROUND PRESENTATION 


A, THE NEED FOR A TUNED MOUNTING FIXTURE FOR THE MWSM 

In an effort to the increase the survivability of ships in 
battle, more attention is being focused on the shock 
qualifications of combat systems equipment installed on ship. 
To ensure the reliability of this equipment in a mechanical 
shock environment it must first be "shock qualified" on the 
appropriate shock machine. For the test to be sufficiently 
valid, the mechanical shock excitation to the equipment should 
be simulated as closely as possible to the expected local ship 
structure response to an underwater explosion. 

MIL-S-901D requires that the test item be mounted to the 
shock machine anvil table in a manner characteristic of its 
shipboard orientation. For most equipment, this requires using 
a specified number of support channels and rails which make up 
the standard fixture (Figure 2). Chalmers and Shaw (1989) note 
that most users of MIL-S-901D believe that using the standard 
fixture on the M('7SM produces a more severe, all encompassing 
shock excitation than would actually be experienced onboard 
ship. They provide evidence in their report that using the 
standard fixture can be an undertest as well as an overtest 
since a high frequency mounting will indeed pass higher 
acceleration levels at higher frequencies, but will not 
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provide the resonant euiplif ications generated by a lower 
frequency fixture. In essence, if the ship's structure excites 
the equipment at a lower frequency than was supplied by the 
MWSM, the installed equipment may experience severe resonant 
vibrations that were not experienced during testing. If a low 
frequency fixture such as the soft deck simulator were used, 


it too 

may be an undertest 

since 

it would 

not excite 

the 

higher 

resonant frequencies 

that 

may be 

present in 

the 


equipment. 

As noted in Corbell(1992), a finite element transient 
shock analysis of the DDG-51 Class Deck House was conducted by 
the Ship Structure and Protection Department of David Taylor 
Research Center (DTRC/UERD). The preliminary report by 
Costanzo and Murray (1991) was obtained along with the 
predicted shock analysis for various weight combat systems 
equipme**: located on the 0-3 level of the DDG-51 class ship. 
From this report, the predicted shock excitations to three 
representative combat systems equipments were used as the 
desired response for a tuned fixture on the MWSM. The 
following equipments selected fall within the weight range of 
test items normally qualified on the MWSM: 

- Radar Receiver Transmitter (RT-1293/SPS-67) 325 lbs 

- Beam Programmer {MX-10873/SPY-1D) 1000 lbs 

- Radio Frequency Amplifier (AM-7159/SPY-1B) 4600 lbs 
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The predicted acceleration wave forms and associated shock 
spectra, analyzed for maximum shock trial severity, for the 
three equipments are show in Figures 4 through 6. 

The shock spectra which are insensitive to small waveform 
variations, describe the characteristic frequencies of the 
shock induced excitations. The shock spectrum is generated by 
plotting the maximum absolute response of a single-DOF 
undamped oscillator as a function of its natural frequency 
when subjected to the base excitation of interest. It is 
important to note, that when attenpting to compare rhock 
motion, the wave "form" in the time domain is much less 
important than the wave “characteristics" displayed in the 
shock spectrum (frequency domain). Therefore, when trying to 
reproduce the shipboard excitation on the MWSM it is not 
necessary to reproduce the wave form, only the frequency 
components at proper acceleration amplitude levels. TVie 
predicted shock spectra demonstrate that the deck house shock 
environment will be significantly different than the high 
energy, single frequency shock excitation produced by using 
the standard fixture on the MWSM. MIL-S-901D does not require 
that any particular waveform or spectrum be reproduced, 
however the need for a timed mounting fixture for the MWSM 
clearly exists. 
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Figure 5. DTRC/UERD Predicted Beam Programmer Vertical 
Acceleration Wave Form and Shock Spectrum. 
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PiSTure 6. uT'-f'/'JERD Predicted Radio Frequency Amplifier 
Vertical Acce ’ e .‘ation Wave Form and Shock Spectrum. 
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B. jmparing shipboard and tuned mounting fixture response 

Extending the concept developed by Chalmers and Shaw 
(1989), Corbell proposed that a two degree-of-freedom (DOF) 
fixture could be designed for the MWSM that would excite an 
attached piece of equipment with most of the energy 
concentrated at two resonant frequencies. If the fixture was 
"tunable", a MIL-S-901D user could refer to UNDEX shoc)c 
spectra data or modal testing data for a specific area of a 
ship and tune the fixture to simulate two of the dominant 
excitation frequencies. 

Figure 7 shows a comparison of the DTRC/UERD DDG-51 
predicted shocV: spectrum for a radar receiver/transmitter to 
the computer modeled shoclc spectra of the Standard Fixture and 
two DOF fixture on the MWSM. The DTRC/UERD predictions show 
significant levels of acceleration amplitudes at 55 Hz and 155 
Hz. The modeled Standard Fixture provides a high level of 
acceleration at 72 Hz but does not supply a sufficient 
amplitude of acceleration between 130 Hz and 170 Hz to 
adequately shoc)c qualify the radar with respect to the 
DTRC/UERD predictions. 

The two DOF model was "tvined" to better simulate the shoc)c 
characteristics predicted for the radar receiver/transmitter. 
As can be seen from Figure 7, the two DOF model reproduced the 
equipment accelerations at the characteristic frequencies and 
cimplitudes of the DTRC/UERD predictions. Equipment qualified 
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Figure 7. Comparison of the DTRC/UEPD Predicted Vertical 
Acceleration Shock Spectra to the Modeled MWSM Fixtures for 
the Radar Receiver/Transmitter - 325 lbs 
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in this manner could be highly relied on to perform in a 
shipboard shock environment. 

C. TWO DEGRE32--OF-FREEDOM FIXTURE MODEL 

1. Coupled and Uncoupled Natural Frequencies 

The Ml'JSM with equipment attached was modeled as two 
mass-spring-danper systems coupled together as shown in Figure 
8. The anvil, when struck by the hammer, experiences a 
half-sine wave vertical acceleration impulse of approximately 
one millisecond in duration. The magnitude of the impulse is 
controlled by the hammer height prior to release. This impulse 
was modeled as shown in Figure 9. 


TWO-DEGREE-OF-FREEDOM FIXTURE .MODEL 



CiLiJ 


C2LiJ 


ANVIL 


K1 



K2 


Ml = Upper Tier Mass 

Includes; Equipment Mass 

Mounting Hardware 
Effective Spring Mass 

M, 3= Lower Tier Mass 
^ Includes: Intermediate Mass 

Effective Spring Mass 

K = Tier Stiffness 
C “ Tier Structtural Damping 




Fisrure 8. Modeled Two Degree of Freedom Fixture Subjected 
to Base Excitation. 


15 


I 









MODELED ANVIL ACCELERATION FOR 1 FOOT HAMMER DROP 



Figure 9. Modeled Anvil Table Acceleration for a One Foot 
Haramer Drop. 


The 

mass Ml 

represents 

the 

equipment mass 

and 

associated 

mounting 

hardware 

while 

Mj represent 

the 


intermediate mass 2 md support movintings. Each tier possesses 
characteristic stiffness and dairying properties designated by 
K cind C respectively. 

To facilitate the solution of the mathematical model, 
the equations of motion for the system can be expressed in 
coordinates relative to the motion of the base, z. Expressing 
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the relative coordinates, and yj in terms of absolute 

coordinates, 

(1) 

y,=Jc^-z 

y2=x2-z 

y-^^x^-t ( 2 ) 

72-^2-i 

yields the relative coordinate transformations for 
displacement, velocity and acceleration respectively. Solving 
for the coupled equations of motion gives the following matrix 
equation: 


k o] 

y] 

* 4* 

O 





-K 

(iCi+iCj) 


0 Mil ( 3 ) 

0 


By assuming the structural damping is small, the 
natural frequencies for free vibration can be found by 
assuming the damped and undamped natural frequencies are 
approximately equal. Neglecting the deimping matrix: 


k o] 


k -^1 

W- 

o' 

o 

k) 

-kj^ {/Cj+iCj) 


0 


( 4 ) 
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The system natural frequencies can then be found by solving 
the following determinate for the system eigenvalues: 


detjf/C] -CillM] ]|=0 


where < 0^=2 


(5) 


Defining the ratio between the upper and lower tier 
masses as. 


a = 


5 : 


( 6 ) 


the relationship between the system's coupled and uncoupled 
natural frequencies^ can be described by, 


V2 

4^= — v'A'f ) "- (2f 1 2 (8) 


^ These equations are based on the valid assumption 
that the structural damping is small and that the damped and 
undamped natural frequencies are equivalent. 
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where f^i and fn 2 are the coupled system natural frequencies 
and and fj are the uncoupled tier natural frequencies 
defined by: 


f =~i_ 

2i^J 




^2 = 


2ll\ 


h. 


(9) 


( 10 ) 


If the desired system natural frequencies are )cnown, 
equations (7) and (8) can be iteratively solved for the 
required uncoupled natural frequencies. Knowing the mass 
associated with each tier, equations (9) and (10) can be used 
to tune the fixture to the required stiffness, K. 

From the shoclc spectrum shown in Figure 7, Corbell 
chose 60 Hz and 155 Hz as the system design resonant 
frequencies of the two DOF model. Figure 10 represents the 
iterative solution of equations (7) and (8) for the uncoupled 
natural frequencies in graphical form. As an example, the data 
was plotted as a function of the mass ratio a=l while fixing 
the lower tier natural frequency at f2=100 Hz. Choosing an 
upper tier natural frequency of 94 Hz from the coordinate axis 
provides the desired system natural frequencies of 60 Hz and 
155 Hz. This procedure was repeated for a mass ratio of 0.7 
and 1.3 as shown in Figures 11 and 12. 
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Figure 10. System Coupled Natural Frequencies as a Function 
of the Upper Tier Natural Frequency (fj . Lower Tier Natural 
Frequency (fj) Fixed at 100 Hz. ioass Ratio = 1. 
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Figure 11. System Coupled Natural Frequencies as a Function 
of the Upper Tier Natural Frequency (f^) . Lower Tier Natural 
Frequency (f2) Fixed at 100 Hz. Mass Ratio = 0.7. 
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Figure 12. System Coupled Natural Frequencies as a Function 
of the Upper Tier Natural Frequency (f;) . Lower Tier Natural 
Frequency(f 2 ) Fixed at 100 Hz. Mass Ratio=1.3. 
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2. Analytical Response Calculations 

To determine the response of the two DOF fixture to 
base excitation, the matrix equation 3 was first uncoupled 
using the normalized transformation matrix, 



( 11 ) 


where the columns of (<t>] are the mode shape vectors 
(eigenvectors) . The resulting uncoupled equations of motion in 
natural coordinates are: 

fii+2CoJfljf|i+o)LTi^ = -2i (12) 

ti2+2CQ^2f|2+wLT)2 = --22 

where ^ is the critical damping ratio estimated from the MWSM 
calibration test data cortpiled by Costanzo and Clements 
(1988). Examination of the calibration data showed that the 
test weight acceleration response damped out in approximately 
.3 to .5 seconds after hammer iirpact. This corresponds to a 
critical dairying ratio of between three to five percent. 
Figures 13 through 15 show the damped acceleration response 
for the two DOF model for three values of the critical danping 
ratio. 


23 



Figure 13. Modeled Two DOF Fixture Acceleration Response 
for Zeta = 0.03. 
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Flgxirs 15. Modeled Two DOF Fixture Acceleration Response for 
Zeta =5 0.05. 
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Equations (12) and (13) where solved independently 
using a continuous-time linear system simulation. The solution 
in natural coordinates was then transformed into relative 
coordinates using equation (11) and then to absolute 
coordinates using equation (1) and (2) to obtain the system 
displacement, velocity and acceleration responses. The MATLAB™ 
algorithm solving the two DOF system is presented in 
Appendix A. 

The two DOF model and numerical response algorithm 
provided the prerequisite information needed for designing and 
predicting the response of the proposed fixture. In Chapter 
IV, a semi-definite three degree of freedom model will be 
introduced which accounts for the unconstrained 4500 lb anvil 
table interaction with the fixture after the time of hammer 
impact. 
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III. TWO DEGREE-OF“FREEDOM TUNED FIXTURE DESIGN 

The designs discussed in this chapter are based on the two 
DOF model presented in Chapter II. Tne procedure followed vas 
to first identify the desired system natural frequencies for 
the fixture and then determine the required tier natur>^x 
frequencies using equations (7) and (8) . Once the tier natural 
frequencies were known, the required stiffness for each tier 
could be determined using equations (9) and (10). The mass, 
stiffness and damping for each tier v;ere input into a 
niimerical algorithm to determine the fixture response to the 
half-sine wave impulse shown in Figure 9, The displacements of 
each tier relative to the anvil table were determined for a 
variety of hammer heights, equipment weights and system 
natural frequencies. The maximum dynamic relative displacement 
within the fixture was obtained for the above conditions and 
then used to conduct a static stress analysis of the fixture. 
A detailed procedure of the design process is presented below. 

A. DESIGN OF THE TWO DOF TUNED FIXTURE 


1. Design Considerations 


The 

dec 

ign 

phase of 
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research began 

with 
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frequencies desired for 
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fixture. After 

referring with 

Mark 

McClean at Naval 

Sea 


Systems Command, the previously desired frequencies of 60 Hz 
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and 155 Hz had been lowered to 30 Hz and 80 Hz for the goal 
frequencies. Modal testing in the applicable areas of the ship 
indicated that the shipboard frequencies were lower than the 
DTRC/UERD predictions. 

A second design consideration was cost. The Navy could 
easily develop a pulse-shape matching machine capable of 
reproducing shipboard shock spectra but the cost per machine 
would be significant. The Navy currently has 13 MWSM in use 
and replacing them with million dollar machines is not 
economically feasible. The design approach was to develop a 
fixture that used many of the current parts and accessories 
already in use on the MWSM. By doing so, fabrication and 
material cost could be minimized. 

MIL-S-901D provides some general guidelines for 
designers building equipment to endure shock excitations, many 
of which are applicable to the fixture design. Materials 
recommended for structural members are those with high yield 
strength, high ductility, high fracture touglmess and when 
possible, light weight. Cast iron and cast aluminum have 
generally proven to be unsatisfactory when used as strength 
members due to their high notch sensitivity cind brittleness. 
Areas normally subjected to conpressive stress under static 
loading conditions may experience tensile stress when 
subjective to the cyclic loading experienced under shock and 
vibration. This fact must be kept in mind specifically when 
designing welded and bolted joints. These joints tend to cause 
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stress concentrations which may be acceptable in compression 
but may lead to failure in tension. 

2. Preliminary Design 

The following designs were based on the two degree of 
freedom model shown in Figure 8. The focus here was to define 
a structure with a suitable set of tier stiffness and mass 
components which would provide the desired response without 
yielding or structural failure. Numerous designs were explored 
and then narrowed to two for this research. 
a. Coil Spring Deck Design 

The soft deck simulator shown in Figure 3 provides 
some obvious quali.ties which would be quite desirable for this 
design. The most significant quality being the ease in tuning 
the fixture. As the weight of the attached equipment varies 
the fixture is easily tuned by adding or removing spring 
cartridges to maintain the desired natural frequency. A second 
quality is the relatively large vertical displacement allowed 
within the fixture. Since the relative displacement is 
inversely proportional to square of the natural frequency, it 
becomes increasingly necessary to allow more relative motion 
between the tiers as lower frequencies are trying to be 
obtained. The coil springs will allow up to two inches of 
travel without being over stressed or becoming coil bound. 

Another feature of using coil springs is the 
ability to distribute the spring force over the entire surface 
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area of the table top which has two advantages. First, as 
mentioned in Clements (15)72), the MWSM is not a perfect 
machine and does not alvays produce a strictly vertical 
impulse on the anvil table. By distributing the spring force, 
the fixture becomes more stable since it can better resist 
the rotational moment produced by the supported mass when its 
center of gravity is not in line with the direction of the 
impulse force. Secondly, since the force is distributed, the 
bending moments applied to the rigid intermediate mass (Mj) 
are minimized. This reduces unwanted secondary frequencies 
generated by deflections within the structure of the 
intermediate mass. 

To achieve the desired system natural frequencies 
of 30 Hz and 30 Hz the required uncoupled frequencies from 
equations (7) and (8) were determined to be fi=48 Hz and f2=50 
Hz. The mass ratio a was set equal to one with each tier mass 
weighing 3000 lbs. This would leave 1400 lbs available of the 
7400 lbs allowable on the MWSM for mounting hardware and the 
inclining fixture called for in MIL-S-901D. Using equations 
(9) and (10) the required upper and lower tier stiffness was 
determined to be Kj=7.66x10® Ib/in and K2=7.06x10® Ib/in. 

The coil springs currently used in the soft deck 
simulator have a stiffness* of 2205 Ib/in and weigh 

* This is the mean stiffness per coil spring as 
determined by Steve Schecter at Hughes Air Craft, Fullerton, 
Ca. 






approximately 15.5 lbs each. This would require 348 springs on 
the upper tier and 320 springs on the lower tier for a 

combined weight of 10,354 lbs for the springs alone. This is 

obviously an unacceptable weight for use on the MWSM. 

An alternate off-the-shelf spring was selected with 
a stiffness of 6090 Ib/in, a weight of 18 lb per spring and 
similar dimensions. This would require 126 springs on the 
upper tier and 116 springs on the lower tier for a total 
weight of 3280 lbs. This reduced the spring weight 
significantly but still consumed too much of the weight 

capacity of the MWSM. The advantages of the coil spring design 
had to be abandoned for a design with a higher 

stiffness-to-weight spring mechanism. The beeun loading concept 
used in the standard fixture was a logical choice. 
b. Beam Spring Design 

■ Obtaining the necessary stiffness while minimizing 
weight encouraged the use of beams for the spring elements. 
The beam elements were modeled as being siimly supported with 
symmetric loading as shown in Figure 16. The stiffness of the 
beam is given by equation (14). 


K=£ = _lil_ 

y a2{3L-4a) 


(14) 
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Simple Supports - Symetrlc Loads 


Figure 16. Beam Bending Model for Determining Beam 
Stiffness . 

This design was based on using the existing 
shipbuilding rails as the spring beam supports. From the 
drawings of the MWSM in MIL~S-901D, this fixed L at 
approximately 50 inches. The load spacing was set to 24 inches 
making a=13 inches. Young's Modules was taken to be E*30xl0‘ 
psi for steel and the area moment of inertia was set to 9.2 
in* , based on the current back-to-back C channels (4x7.25#) 
being used on the standard fixture. The result was a stiffness 
of 2.17x10® lb/in with each beam set weighing approximately 70 
lb. To achieve the desired natural frequencies, the lower tier 
required 4 beam sets ^u:ld the upper tier required 3 beam sets 
for a total weight of 490 lbs. The beams provided a 
significeuit weight savings and became the focus of the 
detailed design. 








B. DETAIL DESIGN OF THE BEAM SPRING FIXTURT^ 

In designing the two DOF fixture, priority was placed on 
the ability to tune the fixture to achieve the desired system 
natural frequencies over the widest possible range of 
equipment weights. An equally important consideration was 
minimizing the weight of the fixture while still maintaining 
the characteristics of the shipboard shock wave. A 1/4 scale 
model of the proposed two DOF fixture is presented in Figure 
17. The f> llowing paragraphs highlight the significant factors 
dealt with during the design process. 

1. Choice of Beams 

The back-to-back C chainnels used or the standard 
fixture were first investigated for use in the two DOF design. 
These channels would mount to the shipbuilding support rails 
in the current fashion using the existing claitping method. 
However, in an effort to simplify the design and to use a 
symmetric cross section, I-beams were ultimately chosen. The 
I-beams used in the design are W4xl3 beams with an area moment 
of inertia of 11.3 in* resulting in a stiffness of 2.46x10* 
lb/in. By modifying the existing clcunping design, the fixture 
could be assembled without requiring emy drilling or welding 
of the I beams in areas of high stress. To achieve the system 
natural frequencies desired, required four I-beams on the 
bottom tier and three on the top tier. To add more flexibility 
in mounting the equipment to the fixture, the design allows 


34 



















for interchanging the I-beams with the back-to-back C channels 
on the upper tier. 

2. Tuning The Fixture 

As the equipment weight attached to the fixture 
changes, the fixture is tuned by altering the 
stiffness-to-mass ratio of the upper euid lower tiers. On the 
upper tier, the stiffness is increased by reducing the 
dimension "a“ shown in Figure 16. The initial design of the 
lower tier allowed for the same adjustment to alter the 
stiffness. In an effort to increase the rigidity of the 
intermediate mass Mj, and to increase stability by lowering 
the center of gravity of the entire fixture, the design was 
changed fixing a=13.5 inches for the lower tier beams. The 
lower tier natural frequency is adjusted by adding or removing 
ballast to the intermediate mass as necessary. 

When tuning the fixture, the effective mass of the 
springs was taken into account. Using Furtis (1972) as a 
reference, the effective mass of the spring elements was 
determined to be 52% of the total beam mass. For K^, this mass 
was added to the intermediate mass Mj. For Kj, 52% of the 
spring mass was added to the upper tier mass and the 
remaining 48% added to Mj since it would be in motion with the 
intermediate mass. The distribution of the effective spring 
mass is some what trivial since its magnitude is quite small 
relative to the tier masses. 




3. Intermediate Mass (M,) Design 

As shown earlier in equations (7) and (8), the ratio 
of the upper tier mass to lower tier mass (a) has an effect on 
the system response to the uncoupled tier frequencies. The 
plot in Figure 18 depicts this relationship. For desired 
system natural frequencies of 30 Hz and 80 Hz, the plot shows 
the required tier natural frequencies as a function of a. As 
a increases, the required upper tier natural frequency 
decreases while the lower tier increases. For a greater than 
1.3, the required uncoupled natural frequencies diverge very 
quickly indicating that the desired system response is not 
achievable. This signifies that the intermediate mass, Mj, 
must be greater than or equal to 77% of the equipment mass. 
The design weight of Mj including the mounting hardware is 
17 00 lbs and will require adding ballast for equipment weights 
over 2200 lbs. The result is, will consume a significant 
portion of the weight capacity of the MWSM and will limit 
equipment weight to approximately 3500 lbs. This weight 
allowance, although somewhat reduced, will still allow 
significant numbers of equipments to be qualified with this 
design. 

A second design consideration for M, was that it must 
be stiff enough to appear as a rigid body to the rest of the 
fixture. This would prevent unwanted deflections and 
vibrations not associated with the desired system response. 
The necessary stiffness was achieved by forming a grid 
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Figure 18. Required Upper and Lower Tier Natioxal Frequencies 
as a Funct:ion o£ the Mass Ratio. 
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Mass Ratio (Mi/M2) 




structure of I-beams with an overall stiffness of l.OxlO’ 
Ib/in which is approximately 13 times higher than Kj or Kj. 

To preclude the need for fabricating a different 
spring beam for the upper tier and lower tier, the 
intermediate mass provides the same mounting geometry as 
the shipbuilding support rails mounted on the MWSM anvil 
table. This allows interchanging the upper and lower tier 
springs as well as the back-to-back C channels currently used 
as part of the standard fixture. In some instances, mounting 
equipment on the upper tier with back-to-back C channels will 
be more convenient and provide more flexibility in cuning the 
fixture. 

C. STRESS ANALYSIS OF THE TWO DOF FIXTURE 

Considering cnly deflections in the vertical direction, 
evaluation of the most likely modes of failure focused on the 
tensile stress in the fasteners and bending stress in the 
spring beams and clamps. The numeric model was loaded to its 
weight capacity and then excited with the maximum expected 
acceleration level. The upper tier weight was set to 3500 lbs 
and the lower tier weight set to 2700 lbs. Using the 300 Hz 
filtered data shown in Figure 19, the maximum expected anvil 
acceleration was picked to be 200 g's for a five foot hammer 
drop. Using these parameters in the numeric model, the 
relative displacements of the upper and lower tiers with 
respect to the anvil table were determined. Figure 20 is a 
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Figure 19. MWSM Peak Anvil Table Acceleration vs. Hammer 
Height. Courtesy of Clemenets (1972). 
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spring beam can be calculated. 


Af=Fa 


(16) 


Knowing the bending moment ,M, the maximum bending stress 
was determined from equation (17); 


o 


t>ending' 


Me 

I 


(17) 


where c=2.0 inches for the four inch I-beam. The shear stress 
in tne beams was determined to be negligible and not included 
in the calculations. (Shear stress is zero at the extreme 
fiber were the bending stress is maximum). 


42 





The equivalent static force ,F, was used correspondingly 
to determine the stresses in the clanps and bolts fastening 
the fixture together. The maximiom stress in the spring beams 
was determined to be 47 kpsi and 13 kpsi in the clamps. 
Maximum tensile stress in the bolts was determined to be 7.5 
kpsi . 

D. DESIGN DRAWINGS AND MATERIAL SELECTION 

Detailed design CAD drawings of the proposed two DOF tuned 
fixture were provided to Naval Sea Systems Command for 
approval and fabrication. A set of these drawings have been 
included in Appendix B of this document. 

The material selected for the spring beams has a yield 
strength of 50 kpsi. Materials selected for all other 
structural elements is standard 36 kpsi yielc* strength mild 
steel. Bolt material is high quality grade 5 carbon steel 
alloy. All aaterials have been received and fabrication 
expected to begin by 1 June 1993. Upon completion, the fixture 
will be sent to Naval Underwater Systems Center, New London, 
Ct. for testing. 

E. PREDICTED RESPONSE FOR THE TWO DOF FIXITTRE 

To predict the performance of the proposed two DOF 
fixture, a sensitivity analysis of the model was performed to 
determine the tunability and the response characteristics of 
the fixture as equipment weight was varied. For the analysis. 
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the desired system natural frequencies were chosen to be 20 to 
40 Hz for fnj 2 ind 7 0 to 90 Hz for f^. The equipment weight 
including the mounting hardware varied from 500 lb to 3500 lb. 

1. Tunability of the Fixture 

To examine the txinability of the fixture, a worksheet 
was produced in Mathcad® which iteratively solves equations 
(7) and (8) for the required uncoupled tier frequencies f^ and 
fj given the desired system natural frequencies f„i amd f ,^2 and 
the mass ratio a. Using the worksheet, the required tuning can 
easily be determined as equipment weight or desired system 
natural frequencies change. A copy of this worksheet is 
provided in Appendix C. Figure 21 shows the achievcible system 
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Figure 21. System Natural Frequencies as a Function of the 
Upper Tier Weight 








TABLE I. TUNING FOR VARIOUS UPPER TIER WEIGHTS. 


Upper Tier 
Weight*- 
(lb) 

Lower Tier 
Weight 
(lb) 

“a" Upper 
Tier^ 

(inches) 

# Beams 
Upper 
Tier 

fni 

(Hz) 

fnj 

(Hz) 

500 

2500 

16 

2 

54 

93 

750 

2500 

16 

2 

49 

88 

1000 

2500 

16 

2 

45 

83 

1250 

2500 

16 

2 

42 

80 

1500 

2500 

16 

2 

38 

79 

1750 

2500 

16 

2 

36 

78 

2000 

2500 

16 

2 

34 

77 

2250 

2500 

16 

2 

32 

76 

2500 

2500 

15 

2 

32 

77 

2750 

2500 

15 

2 

30 

77 

3000 

2700 

14 

3 

33 

84 

3250 

2700 

14 

3 

32 

83 

3500 

2700 

14 

3 

31 

83 


-Upper tier weigiic includes equipment, mounting hardware and 
effective spring weight. 

2-Upper tier equiment mounting spacing. Refer to ficpare 16. 


natural frequencies as the weight of the upper tier is varied 
between 500 and 3 500 lbs. TABLE I shows the necessary 
adjustments which would be required to tune the fixture. The 
sharp rise in system natural frequencies for upper tier 
weights less tham 1500 lbs indicates that adding ballast to 
the upper tier may be required for some light weight 
equipment. 
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A second analysis was performed to determine the range of 
system natural frequencies obtaineOole for a specified 
equipment weight. The equipment weight was given to be 2500 
lbs and the fixture tuned by varying the weight of Mj and 
adjusting the equipment mounting spacing ("a" in Figure 16). 
For an equipment weight of 2500 lbs. Figure 22 demonstrates 
the obtaincible system natural frequencies. The data provided 
in TABLE II indicates the necessary adjustments. Figures 21 
and 22 demonstrate that this fixture design could be tuned to 
the desired system natural frequencies under a variety of 
loading conditions. 
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TABLE II. TUNING FOR AN UPPER TIER WEIGHT OF 2500 LBS. 


Case 

Number 

Lower Tier 
Weight 
(lb) 

"a" Upper 
Tier 
(inches) 

# Beams 
Upper 
Tier 

fn, 

(Hz) 

fnj 

(HZ) 

1 

3000 

16 

2 

31 

70 

2 

3000 

15 

2 

31 

71 

3 

2750 

16 

2 

31 

73 

4 

2750 

15 

2 

32 

75 

5 

2500 

16 

2 

31 

76 

6 

2250 

15 

2 

32 

81 

7 

2250 

14 

2 

33 

83 

8 

2000 

13 

2 

34 

89 

9 

2000 

14 

3 

37 

96 

10 

1750 


3 

38 

107 
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2. Two DOF Fixture Response Characteristics 

As mentioned in Chapter II, for the fixture to be 
effective, it should provide the equipment with accelerations 
that are "characteristic" of the expected shipboard 
excitations. Since the desired frequencies have been lowered 
by NSWC to 3 0 Hz cind 80 Hz, it is difficult to compare the 
modeled results to the DTRC/UERD predictions. However, as 
shown for five different equipment weights. Figures 23 through 
2 8 demonstrate that the modeled time history accelerations and 
associated shock spectra provide significant amplitudes of 
acceleration at the desired frequencies. 
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ACCEI.EKATION (g’s) ACCELERATION (g's) 


MODELED TWO DOF FDCTURE UPPER TIER ACCELERATION 



modeled TWO DOF FDCTURE ACCELERATION SHOCK SPECTRA 

450,----- 



Figure 23. Predicted Acceleration Wave Fonti and Shock 


Spectrum for 1500 lb Upper Tier Weight. 
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ACCELERATION (g‘s) ACCELERATION (g’s) 


MODELED TWO DOF FDCTURE UPPER TIER ACCELERATION 



MODELED TWO DOF FIXTURE ACCELERATION SHOCK SPECTRA 



FREQUENCY (Hz) 


Figure 24. Predicted Acceleration Wave Form and Shock 
Spectmm for 2000 lb Upper Tier Weight. 
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ACCELERATION (g's) ACCEI.ERATION (g's) 



Spectrtim for 2500 lb Upper Tier Weight. 
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ACCELERATION (g's) ACCELERATION (g’s) 



Spectrum for 3000 lb Upper Tier Weight. 
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ACX:ELERATJ0N (g’s) ACCEI.ERATJON (g's) 


MODELED TWO DOF FIXTURE UPPER TIER ACCELERATION 

5 Foot Hammer Drop (200 g's) 

_ Upper Tier Weight = 3500 Ib 
Lo^r Tier Weight = 2700 lb 
Zeta = 0.04 
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Figure 27. Predicted Acceleration Wave Form and Shock 
Spectrum for 3500 lb Upper Tier Weight. 
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IV. SEMI-DEFINITE THREE DEGREE-OF-FREEDOM MODEL 

As mentioned in Chapter II, Corbell formulated a two- 
degree-of-freedom model to simulate the fixture response to 
the anvil table input. The fixture model was shown in Figure 
8. The input into the model was the half-sine wave base 
acceleration impulse described in Figure 9. After the base 
impulse passed {1 msec), the anvil table was considered 
stationary and the fixture vibrated as a two DOF mass-spring- 
demiper system attached to a fixed foundation. This model 
provided useful information required for determining the 
feasibility and design of the proposed fixture. In this 
chapter, a semi-definite three DOF model will be introduced 
which will take into account, the 4500 lb anvil table 
interaction with the fixture after the time of hammer impulse. 

A. SEMI-DEFINITE THREE DOF MODEL DERIVATION 
1. Model Description 

As presented in Clements (1972), the anvil table is 
bolted to the machine foundation with 12, 2-inch-diameter 
bolts in a manner which allows the table three inches of 
vertical travel after hammer impact. This vertical distance 
can be decreased by raising the table with pneumatic jacks 
prior to performing the test. Approximately 50 msec after 
hammer impact, the table reaches the limit of vertical travel 
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and experiences a sharp "table reversal" in the form of a 
fairly simple negative acceleration pulse. Since the duration 
of the reversal is 2 to 4 times that of the hammer impact and 
energy has been expended ir the system, the magnitude of the 
acceleration pulse is much less than that caused by the 
hammer. As the table falls back down on the foundation it 
experiences another impulse of longer duration and 
consequently lower acceleration. From the time of hammer 
impact to the time of table reversal, the anvil table is 
unconstrained and the entire structure behaves as the semi- 
definite three DOF system shown in Figure 28. Considering the 
initial impact as significantly more sever than the 
subsequent, it was possible to simplify the model and not 
include the table stop in the mathematical analysis for 
determining the characteristics of the fixture response. 

Referring to Figure 28, the "upper tier" mass Mj 
represents the mass of the equipment and mounting hardware and 
the effective spring mass for the upper tier. represents 
the intermediate mass and the effective spring mass for the 
lower tier. Mj is the mass of the anvil table. The spring 
stiffness is designated by Kj for the upper tier and Kj for the 
lower. Again, as in the two DOF model, the structural damping 
is considered proportional to the mass and stiffness. The 
upper tier and lower tier damping is designated by Cj and Cj 
respective .'.y. 
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THREE DOF MODEL 



Three DOF Model for Fixture on MWSM 








Referring to Clements (1972), the peak anvil table 
acceleration produced by the half-sine impulse is very nearly 
a linear function of the hammer impact velocity and is 
essentially "independent" of the load when it is channel 
mounted. Using this information, the forcing function F(t) can 
be modeled as simply the anvil acceleration times the anvil 
mass and is independent of the load attached. 


Fit) (18) 

Since the test scnedule for the MW3M listed in MIL-S- 
901D specifies the "hammer drop height" for a given table 
weight, it is convenient to write F(t) as a function of the 
height of the hammer drop. Referring back to Figure 19, the 
300 CPS filtered data shows a linear correlation between anvil 
acceleration and velocity with a slope of approximately 11.35 
g'sec/ft. Using the relation, V=V2gh, the peak anvil table 
acceleration can be related to the hammer height by the 
followii relation: 


Ji^ipeak) =(11.35) (32.2)/2^ 
where: h=Haininer Height Above Anvil Table in Feet 


(19) 


57 




Knowing the peak acceleration of the half sine wave 


impulse, the forcing function F(t) can then be modeled as; 


F{C] =293 3v^(Ar3) sin (G)t) 0;Jti0.001 sec 

F(t) -0 t>0.001 sec (20) 

where: <.)=l000*r 

The 1000 CPS filtered data shown in Figure 19 
represents the table’s 750 Hz longitudinal mode of damped 
vibration which last for only about 5 periods and does not 
significantly contribute to the motion of the fixture. 

2. Mathematical Analysis of the Three DOF Model 

The number of degrees of freedom of a system is 
defined as the number of independent coordinates necessary to 
describe the motion of a system completely. For the model 
shown in Figure 28, this will recruire three independent 
coordinates to define the motion of the three masses of the 
system. These coordinates, Xj, and Xj, define the absolute 
displacements of the masses, Mj, Mj and Mj, respectively. The 
time derivatives of these coordinates yield the velocity and 
acceleration: 


1^2 ' ^2 

(Displacement) 


,X2 

(Veloci ty) 

( 21 ) 

1 ^2 1 ^3 

(Accel elation) 
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The resulting equations of motion in absolute 
coordinates can be written as: 


^kj^x^-c^X2 =0 (22) 

1112X2 + (Cl + c- ) X2 + ( +;:2 ) -^2 “ c^iXi - JCjXi - C2X3 -/C2X3 = 0 ( 23 ) 

^ 3 X 3 + 03 X 3 *k2X2-C 2 X 2 -k2X2 =F (24 } 

The above equations of motion can be conveniently 
expressed in the following matrix form: 


■^1 

0 

O' 






O' 





O ' 



0 ] 

0 


0 


^2 

. + 

-^1 

(C1+C2) 

-<^2 

k 

. -f 

-K 

{ k ,^ k 2 ) 

-^2 

^2 


0 

0 

0 

/TJ3 


^ 3 , 


0 

-^2 

•^2 

P 3 


0 

-^2 

K . 

P 3 




Since the structural damping of the fixture is small, 
the system natural frequencies (free vibration) can be found 
by assuming that the damping is negligible and that the damped 
and undamped frequencies are equal. Neglecting danping and 
assximing free vibration, equation (25) can be rewritten as: 


Q O' 



k^ -k^ 0 



o' 

0 

0 

*2 


-ici ( /cj +ic2) -k2 

X 2 

= 

0 

^ 0 0 ^3 

.^3, 


0 -iCj k2 

;^3 


0. 


(26) 
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The system natural frequencies are then found by 
solving the following determinant for the system eigenvalues: 


det|[/0 ]| = 0 

where (j3„=2nf„ 


(27) 


The solution to equation (27) is shown below in 
equations (28) through (30) revealing that one of the three 
system natural frequencies is zero. This zero frequency 
corresponds to che rigid body mode of vibration. 


{Rigid Body Mode) 


(28) 


0 } 


n2' 


1 

B-* 

^| 2 

N 






(29) 


^n3~ 




1 

B-*-. 

2 

N 








( 30 ) 


where 



M," M2 
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As described in Chapter II, if the desired system 
coupled natural frequencies fnj.f^j and f^j are known, equations 
(29) and (30) can be iteratively solved to determine the 
required uncoupled tier natural frequencies fi and f 2 . A 
Mathcad® program is listed in Appendix C which, given the 
c’-'sired system natural frequencies v/ill solve for the required 
tier frequencies . 

Having found the system natural frequencies, the 
natural modes of vibration (eigenvectors) can be solved for 
using equation (31). 

-^1 C 

0 -iCj 

i=l,2,3 

Equation (25) can be transformed into natural 
coordinates by introducing the transformation matrix [<I>] where 
the coluitvns of [<P] are the eigenvectors found from equation 
(31). Using the transforming matrix; 


4>2,i 


= 0 


(31) 


’fl 


4'll 4 ’i2 4>i3' 




= 

4'21 4>22 <l>23 

^2 

(32) 

^3 


4»31 4>32 4>33. 

^3, 
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equation (25) can be decoupled by taking advantage of the 
orthogonality of the modal vectors. Applying the 
transformation matrix and premultiplying by equation 
(25) can be rewritten as: 




Assuming that the modal damping is proportional to the 
mass and stiffness, equation (33) can be decoupled and written 
as three independent differential equations of motion as shown 
below: 


fii+2Co)„iiii+Qfljili=Fi 

(34) 

fi2+2C0)„2f|2+0)L'n2=-F’2 

(35) 

fi3+2{0)„3f|3+W^3Tl3=F3 

(36) 


Equations (34) through (36) were solved for 
displacement, velocity and acceleration using the MATLAB™ 
program provided in Appendix D. As in Chapter II, the critical 
damping ratio was assumed to be 0.04. The solutions in natural 
coordiiiates was transformed back to real coordinates by 
reapplying equation (32). 
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B. THREE DOF FIXTURE MODEL RESPONSE 

The upper tier acceleration response to the half sine wave 
impulse was modeled for the same weight and stiffness 
configurations as presented for the two DOF model in Figures 
23 through 27 . The modeled acceleration time response and 
shock spectra for five different upper tier weights are shown 
in Figures 29 through 33 . As can be seen from the figures, the 
system natural frequencies for the three DOF model are 
approximately 15 to 20 Hz higher than for the two DOF model 
and peak accelerations are less. Figures 34 through 36 show 
the time response and shock spectra comparisons of the two DOF 
and three DOF models. 

Since the anvil table is not absolutely unconstrained due 
to friction on the bolts and contact with the stops, it is 
expected that the actual system natural frequencies will fall 
somewhere between those of the two DOF and three DOF models. 
In the next chapter, the results of the 1/4 scale model 
testing will be presented which will help better predict the 
response of the actual fixture. 
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ACCELERATION (g’s) ACCELERATION (g’s) 


40 


MODELED THREE-DOF FDCTURE UPPER HER ACCELERATION 



Figure 29. Three DOF Model Acceleration Wave Fonti and Shock 
Spectrum for 1500 Lb Upper Tier Weight. 
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ACCELERATION (g’s) ACCELERATION (g’s) 



Figure 30. Three DOF Model Acceleration Wave Form and Shock 
Spectrum for 2000 Lb Upper Tier Weight. 


65 











ACCEl.ERATION (g’s) ACCELERATION (g’s) 


MODELED THREE-DOF FDCTURE UPPER TIER ACCELERATION 


5 Foot Hammer Drop (200 g’s) 
Zeta = 0104 


Upper Tier Weight = 2500 lb 
Loivcr-Tier-Weight- “-2500 Ib- 
Anvil Weight = 4500 lb 


TIME (sec) 


6 C.18 


MODELED THREE DOF FIXTURE ACCELERATION SHOCK SPECTRA 


40 60 SO 1 

.00 

FREQUENCY (Hz) 



Figure 31. Three DOF Model Acceleration Wave Form and Shock 
Spectrum for 2500 Lb Upper Tier Weight. 
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ACX'EI.ERATION (g’s) ACCEl-ERATION (g’s) 


30 


MODELED rriREE-lX)F FDCTURE UPPER TIER ACCELERATION 



Spectr\im for 3 000 Lb Upper Tier Weight. 
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Figure 33. Three DOF Model Acceleration Wave Form and Shock 
Spectrum for 3500 Lb Upper Tier Weight. 
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ACCEl£RAT10N (g’s) ACCELERATION (g’s) 


COMPARING 2 DOF AND 3 DOF MODEL TIME RESPONSE 



TIME (See) 


COMPARING 2 DOF AND 3 DOF MODEL SHOCK SPECTRA 



FREQUENCY (Hz) 

Figure 34 . 2 DOF and 3 DOF Model Time Response and Shock 


Spectra Comparisons for 2500 Lb Upper Tier Weight. 
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ACCELERATION (g’s) ACCELERATION (g’s) 


COMPARING 2 EXDF AND 3 DOF MODEL TIME RESPONSE 


5 Foot Hammer Drop (200 g’s) 

- 2xt3.--iy.04 - 

Upper Tier Weight = 3000 lb 
Low er Tit-r Weirjit = 2700 lb 


\ M 


— 2 DOF 

— 3 DOF 


; ^ f 


j__!_ 


:1 L: 1-' U; L - I 

• i , J/ I 


0 0.05 0.1 0.15 0.2 a25 0_5 0J5 0.4 0.45 

TIME (Sec) 


COMPARING 2 DOF ANT) 3 DOF MODEL SHOCK SPECTRA 

5 Foot Hammer Drop (200 g’s) 

- 1 -- 1 -Upper-Tier-Weigin- =-3000- lb- 

;; 2 DOF Lower Tier Weight = 2700 lb 


-.-r -*-3 DOI^ 

Ml-. 


FREQLTENCY (Hz) 


Figure 35. 2 DOF and .3 DOF Model Time Response and Shock 

SpecLra Comparison for 3000 Lb Upper Tier Weight. 
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FigT.ire 36. 2 DOF and 3 DOF Model Time ReGpo.nse and Shock. 

Spectra Compari.son for 3500 LL Upper Tier VJeight . 
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V. 1/4 SCALE MODEL TESTING AND RESULTS 


To verify the mathematical simulation, a 1/4 scale model 
of the fixture was fabricated for shock testing. The principal 
intent of the testing was to confirm the analytical methods 
employed were appropriate for the application and that the 
expected frequencies characteristics could be obtained. The 
following sections describe the test procedure used and the 
results obtained. 

A. EXPERIMENTAL SET UP 

To supply the inpulse excitation to the model, a drop 
table was constructed which allowed dropping the fixture model 
onto a load cell from various heights. The range of the load 
cell was 0-5000 lbs with a sensitivity of 0.106 mv/lb. The 
upper tier accelerometer had a range of 0-500 g's with a 
sensitivity of 9.98 mv/g. Using a Hewlat Packard 3562A D^Tiamic 
Signal Analyzer, the transfer function between the input 
excitation and upper tier acceleration was obtained. The 
experimental set up is shown in figure 37 and a schematic of 
the instrumentation is shown in figure 38. 

B. 1/4 SCALE MODEL 

The fixture model was built as clo.sely as possible to a 
1/4 scale version of the proposed fixture. Ebcceptions 





































included using aluminum instead of steel for the spring I- 
beams and an additional 1/4 scale reduction in the weight of 
the tiers. The upper tier weight as tested was 43.5 lbs, the 
lower tier weight was 55.2 lbs and the anvil table was 109.1 
lbs. The per-beam-stif'ness of the aluminum I-beams was 
2.09x10* Ibs/in. 

To help verify the boxindary conditions used in the 
mathematical model, the spring beams were mounted using the 
scime clamping design proposed for the full scale prototype. 
The following exaimple provides the dimensions and calculations 
for determining the expected system natural frequencies for 
the three DOF model. 


Example: Detemiine System Natural Frequencies For 1/4 Scale Model Fixture: 
A. SPRING STIFFNESS: 

The spring stiffness is based on the following beam model; 

y 

A 

Simple Supports - Symetric Loads 



1) Dimensions: 

Length: 

Load Applicaton Point; 
Area Moment of Inertia: 
Young's Modulus: 


Upper Tier 

L j =12.5 in 
a I = 3.25 ui 
I I =0.045 in* 
E = 10 10* psi 


Lower Tier 

L 2 =12.5 in 
et =3.25 in 
1 2 - 0 045 in* 
E =10 10* ps; 
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2) Stiffness: 


From the beam bending equation: 




Let x=a and F=W/2; 


W'fl 

yAB*-^ (4 a- 3 L) 

^ 12EI 


Solving for the stiffness: 


, 12 E I ]b 

yAB a^(3L-4a) ^ 


For: 


Upper Tier 


Lower Tier 


(Convert inches to feet) KB j 


12E I 


a j i 3 ■ L I - 4 a j i 


12 KB 


12 E K 


a 2 ; 3 L 2 - 4 a 


(per beam) 

K3 ■ =2.504-10*- 
tt 

KB, =2.504-10^- 
ft 


Upper Tier 

Lower Tier 

Input Number of Beams Per Tier: Beams j = 2 

Beams 2 = ^ 


K J = KB J Beams • 

K 2 = KB 2 Beams-) 

Tier Stiffness: K,=5.008-lo’ - 

' ft 

K, =5.008-10’ - 

ft 

B. SYSTEM MASSES; (Includes Effective Sprina Mass) 


Upper Tier Weiaht; (W'^) 

Intermediate Weiaht: (Wj) 

Anvil Table Weiaht: IV7-^) 

W, =43.5 

W 2 =55.2 

W 3 = 109 1 

M , =—!- 
‘ 32.3 

W 2 
=—=- 

^ 32.2 

W 3 

M 1 = — 

32.2 

M, =1.351 lb—- 
‘ .2 

M 2 ==1714 Ib-^ 

M 3 = 3 388 lb 


sec 


sec 


ft 

sec‘ 
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C. TIER NATURAL FREQUENCIES: 


D. System Natural Frequencies: 

Defining B as; 


Upper Tier 


F =J_. 

F , » 96.904 Hz 


Lower Tier 


“ 2 rzjM 2 


F 2 = 86.023 Hz 


B =■ 


K, K, K,-K, 


M 1 M 3 


M., 


, i| 1 I , 4.K,K, I 

f . =J_. il|B- !B‘-1 —i_.;m , -M-,3', I 

ni I 1 V/ v.< v.< \ 1 - ■’/ i 


N 


M [ M 2 M 3 


First Frequency; 


f^l =77.762 Hz 


' 02 


I 'll '2 2 

‘ ■' -B- IB^- ‘ 


^ M 1 M 2 M 3 


. (M 1 - M 2 “ M 3^, i 


Second Frequency; fj^ = 14^ y 44 Hz 


C. EXPERIMENTAL RESULTS 

The first experiment performed was determining the static 
stiffness of an assembled tier containing three spring beam 
elements. Using the simply supported beam model, the expected 
stiffness of the tier was calculated. The tier was then place 
on a load machine and loaded up to lOOC lbs while measuring 
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the static deflection. Figure 39 shows the analytical 
stiffness versus the experimental. From the plot, the actual 
stiffness was measured to be 43,500 Ib/in which is about 5 
percent higher than the analytical stiffness. This result 
tends to justify the simply supported boundary conditions used 
in the numerical model. 


EXPERIMENTAL VS. CALCULATED TIER STIFFNESS 

0.02S 

0.02 

0.013 

Displacement 
(Inches) 

0.01 

0.005 

0 

200 300 400 500 600 700 *00 900 lOOO 1100 

LOAD (Lbs) 

Figxtre 39. Experimental Tier Stiffness versus Analytical. 

The model was then installed on the drop table to 
determine its response to a six inch table drop. Using the 
load cell as the input and upper tier acceleration as the 
output, the tremsfer function for a ten drop average was 
computed using the HP 3562A. The acceleration data was 
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filtered using an exponential window with a period T=400 msec. 
The data collection was triggered off the load cell input 
using a -5 msec trigger delay. The transfer function with 
phase shift is shown in Figure 40. The experimental transfer 
function reveals distinct peaks at 80 Hz and 160 Hz with 
associated 180 degree phase shifts. These results agree very 
closely with the mathematical predictions of fni=77.8 Hz and 
f „2 = 148 Hz. 

Figure 41 shows the \infiltered load cell and acceleracion 
output sampled at 10,000 Hz. The load cell data shows a 
maximum impulse of approximately 35,000 lbs with a duration of 
about 10 msec. The corresponding upper tier acceleration is 
17 0 g's. Usi.ng a logarithmic decrement as an approximation, 
the critical damping ratio was determined to be two percent. 

The experimental results, although not extensive, support 
the mathematical model very well. They demonstrate that using 
the proposed fixture will transform the high frequency, high 
impact energy of the hcimmer-anvil impact into discrete 
equipment excitation frequencies, tunable to shipboard 
conditions. With full seeVe fixture testii g beginning in July 
1993, promising results are exjjected. 
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Accelerometer Output. Sample Frequency 10,000 Hz. 
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VI 


CONCLUSIONS AND RECOMMENDATIONS 


The mathematical model and experimental results provided 
in this study demonstrate great promise in the concept of 
using a tuned two DOF equipment mounting fixture on the MWSM. 
Using frequencies provided by NSWC, the mathematical model 
provided the required stiffness and mass parameters necessary 
to design the fixture. The design was then scaled to one 
quarter and an experimental model constructed for impact 
testing. The results of the model testing agreed extremely 
well with the results predicted from the mathematical 
modeling. The following is a list of specific comments 
concerning the tuned fixture: 

- The fixture designed in this study is not expected to be 
an all inclusive tuned fixture capable of testing all 
mediumweight equipment. Its specific design is to provide 
excitation at 30 Hz and 80 Hz for equipment up to 3500 
lbs. As more information becomes available on localized 
ship structure excitation, a series of tuned fixtures can 
be constructed (economically), each designed for a 
specific location on the ship. 

- Since the weight capacity of the MT-JSM is limited to 7 4 00 
lbs, equipment weight using a two DOF fixture will be 
limited to approximately 3500 lbs. The ratio between the 
upper and lower tier mass is limited requiring them to be 
nearly equivalent in magnitude. From the DTRC/UERD 
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predictions, this should not be a significant problem 
since heavier equipment will not generally require a two 
DOF fixture to simulate shipboard shock. 

- The I-beam spring design will work very well for 
frequencies above 25 Hz, but would not be suitable for 
lower excitation frequencies. Since the spring beams are 
limited to approximately sixty inches in length on the 
MVJSM, the vertical tier displacements required for lower 
frequencies would result in excessive bending stress in 
the I-beams. For this reason, it is recommended that the 
coil spring design be reinvestigated tor lower 
frequencies. 

- On completion of the prototype fixture in June 1993, a 
series of calibration tests will be performed to determine 
the actual response of the fixture. Based on these test, 
the mathematical model can be updated with structural 
dcimping characteristics as well as adjusting the boundary 
conditions if necessary. This new model will then be 
available for the design of future fixtures. 

The use of the two DOF equipment mounting fixture for 
shock qualifying equipment on the MWSM is in line with the 
Navy's goal of shock hardening ships. Incorporating this 
concept into the current shock qualification procedures w’ill 
provide increased reliability of essential equipment subjected 
to shock from underwater explosions. 
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APPENDIX A. TWO DOF TUNED FIXTURE PROGRAM 


% Program TWODOF.M 

% Lt. David M. Cox 

% Naval Postgraduate School 

% This MATLAB program plots the relative displacement, 
%relative velocity and absolute acceleration of a TDOF mass, 
%spring, damper system subjected to base excitation (half-sine 
%pulse). Additionally, the shock spectra of an undamped SDOF 
%system is calculated using the TDOF upper tier acceleration 
%as the base excitation to the SDOF system. 

clear 


%User Defined Variables 


G=203 

Wtl=3500 

Wt2=2700 

kb=2.6e6 

kl=3‘kb; 

k2=4*kb; 

2 eta=.04 


% Acceleration Magnitude in g's 
% Weight of Upper Tier in Lbs 
% Weight of Lower Tier in Lbs 
% Stiffness per Beam 
% Upper Tier Stiffness in Lbsft 
% Lower Tier Stiffness in Lbs/ft 
% Critical Damping coefficient 


% 

% 

% 

% 


%Formu].ate the Base 
dt=0.0001 
w=1000*pi 
N=5000; 

T=N*dt 

t=[0;dt:(N-1)*dt] 
for i=l:N 
if t(i)<=0.001 

anvil(i) = G*32.2*sin(w*t 
else 

anvil(i) = 0 ; 

end 

end 


Acceleration 

% time step for sampling 
base excitation frequency 
Number of steps 
Total time record 
Time vector from t=0 to t=T 


% Loop to produce half-sine wave pulse 


i) ) 


%Calculate the mass and mass ratio 
ml=Wtl/32.2 ; 
m2=Wt2/32.2 ; 
m=[ml,0;0,m2]; 
alpha=ml/m2; 
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%Calculate Circular Tier Natural Frequencies 
wl=sqrt(kl/ml); 
v2-so''t (k2/m2) ; 

%Systein Natural Frequencies 

var=sqrt ( (v;1^2 + alpha*wl^2 + w2''2)^2 - (2*wl*w2 ) ^^2) 
wnl= (1/sqrt (2 ) ) *sqrt (wl"^2 + alpha*wl^2 + w2^2 - var) 
wn2= (1/sqrt (2 ) ) *sqrt (wl''2 alpha*wl^2 + w2''2 + var) 
%Uncouple the equations 
ull=l; 

u21=(kl - wnl''2 * ml)/kl; 
ul2=l; 

u22=(kl - wn2^2 * ml)/kl; 

U- [ull,ul2;u21,u22] ; 

%Decouple the system 
M=U'*m*U 
%£tiffness Matrix 

k=[kl,-kl;-kl,kl+k2]; % coupled stiffness 

K=U'*k*U 

%Force Coefficients 

FC=(U'*[-ml;-m2]) 

%Solve the uncoupled Equations Of Motion 

% {nddotl} + [2 *zeta*Vi'nl ] * {ndotl} + [wnl^2 ] * {nl} =FC11*anvil /Ml 
% {nddot2) + [2 *zeta*wn2]*{ndotl}+[wn2 ^2]*{n2}=FC21*anvil/M2 

%State Space Matrix for Equation #1 

A1=[0,1; - (wnl'^2) , -2*zeta*wnl] 

B1=[0;FC(1;1)/M(l:l)]; 

Cl=tl,0;0,l;- (wnl''2) , -2*zeta*vml] ; 
Dl=t0;0;FC(l:l)/M(l:l)]; 

[yl]=lsim(Al,Bl,Cl,Dl,anvil,t); 

%State Space Matrix for Equation #2 

A2 = [ 0,1; - {wn2''2 ) , -2*zeta*wn2 ] 
B2=[0;FC(2:2)/M(4:4)]; 

C2 = [ 1,0 ; 0,1; - (wn2''2) , - 2*zeta*vm2 ] ; 

D2=[0;0;FC(2:2)/M(4:4)]; 

[y2]=lsim(A2,B2,C2,D2,anvil,t); 

%Couple The Equations 

Y1DIS=(ull*yl{;,1) + ul2*y2(;,1))*12; 

YlVEL=ull*yl(:,2) + ul2*y2{:,2); 

YlACC=ull*yl(:,3) + ul2*y2{;,3); 

%Absolute acceleration in g's 

X1ACC=(YlACC+anvil')/32.2; 

%Convert to in. 

Y2DIS={u21*yl(:,1) + u22*y2(:,l))*12; 

Y2VEL=u21*yl(:,2) + u22*y2(:,2); 

Y2ACC=u21*yl{:,3) + u22*y2(:,3); 

X2ACC={Y2ACC+anvil’)/32.2; 
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plot(t,XIACC),grid 

title('MODELED TWO DOF FIXTURE UPPER TIER ACCELERATION') 

xlabeK'TIME (sec)') 

ylabel('ACCELERATION (g''s)') 

gtext('5 Foot Hammer Drop (200 g''s)') 

gtext('Upper Tier Weight = 3500 lb') 

gtext('Lower Tier Weight = 2700 lb') 

gtext('Zeta = 0.04') 

meta fig28.met 

plot(t,YIDIS,t,Y2DIS).grid 
xlabel('TIME{Seconds)') 

ylabel ('RELATIVE DISPLACEMFAT' (Inches)') 

title('MODELED TIER DISPLACEMENT RELATIVE TO THE ANVIL TABLE') 

gtext('Hammer Height = 5 feet (200 g'‘s)') 

gtext('1 - Upper Tier Weight = 3500 lb') 

gtext('2 - Lower Tier Weight = 2700 lb') 

gtext('Zeta = 0.04') 

gtext(' 1 '),gtext(' 2 ') 

meta fig 20 .met 

pause 


%Calculate the shock spectr 
NF=120; 

DF=1; 

SF=1 ; 

Freq=(SF;DF:NF]; 
wn=2*pi*Preq; 

F=[0;-1]; 

GG=tl,0]; 

H=l-1] ; 


% Number of Frequencies 
% Frequency Increment 
% Start Frequency 
% Frequency Vector 

% State Space Matrices 


for i=l:NF 


% Step Thru Natural Frequencies 


E=[0,l;-(wn(i)'2),0]; 

GG=[wn(i) ''2,0]; 

[yspec]=lsim(E,F,GG,H,XIACC,t); 
xspec = yspec + XIACC; 
maxspec(i)=max(abs(xspec)) ; 
end 


plot(Freq,maxspec) 

title('MODELED TWO DOF FIXTURE ACCELERATION SHOCK SPECTRA') 
xlabel('FREQUENCY (Hz)') 
ylabel('ACCELERATION (g''s)') 
grid 

gtext('fnl = 38 Hz') 
gtext('fn2 = 78 Hz') 
meta fig28.met 
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APPENDIX B. TUNED FIXTURE DESIGN DRAWINGS 


DESIGN DRAWINGS 


ITEMIZED PART LIST 


1 Piece Number 

Item 

Size (inches) 

No. Required 

' 1 

Shipbuilding Channel 

Ixll.lU 

2 

'l 2 

I-Beam 

W5 X 18.5# 

4 

i 3 

I-Beam 

W8 X 67# 

3 

1 ^ 

Block 

l"x2 1/4" X 10" 

28 

1 5 

Stiffner 

See Figure 2 Sheet 2 

6 

6 

I-Beam 

W4x 13# 

7 

1 - ^ - 

Stiffner 

See Figure 2 Sheet 5 

12 

1 8 

Block 

I"x2"x5 1/8" 

28 

1 9 

Clamp 

See Figure 4 Sheet 1 

14 

1 10 

Clamp 

See Figure 5 Sheet 1 

28 

i_y_ 

Washer 

1 l/8"IDx2"OD 

84 

^ m. 

Hexagon Head Bolt 

l"-8x4 1/4" Long 

84 


Notes: 

1) Top flange of piece nurnber 1 shall be burned off or 
cut off to a width of 1 3/4". 

2) Piece numbers 9 and 10 will be shaped to fit inner 
surface of piece numbers 1 and 6 respectively. 

0 ) Hexagon head bolts may be replaced with eq\iivalent 
alien head bolts ease installation of piece 10 to 
piece 2. Depending on weld thickness and assembly 
tolerance, use of a socket to tighten hexagon head 
bolts may be difficult. 
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APPENDIX C. MATHCAD® WORKSHEETS 


APPENDIX C1 

MathCad Worksheet To Find Two DOF Fixture Ucoupled Frequencies: 


input the Desired Coupled Natural Frequencies: (Hz) 

=30.00 --80.00 

Input The Mass Ratio: 

Mass Ratio: a = I 

Initial Guess at uncoupled frequencies: f ] =40 f 2 =40 

Given the initial guess for and f2, the following equations are iteratively solved for the required 
values of fi and f2. 

Given 



SOLUTION 

Required Uncoupled Tier Natural Frequencies: f j = 49.%3 f 2 = 48 036 Hz 
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APPENDIX C2 


Mathcad Worksheet to Determine System Natural Frequencies Of Two DOF Fixture: 


A. SPRING STIFFNESS: 


The spring stiffness is based on the following beam model: 

y 



Supply the following Dimensions; 


1) Dimensions: 

UDoerTier 

Lower Tier 

Length: 

L ] = so in 

L 2 = 50 in 

Load Application Point: 

a I = 13 in 

8 2=13 in 

Area Moment of Inertia; 

1, ^113 in* 

I 2 =113 m* 

Young’s Modulus: 

E =30 10* psi 

E =30 10* psi 


2) Stiffness: 


From the beam bending equation; v (x^ + 3 - 3 l) 

' ^ 6 El 


Let x=a and F=W/2; 


W 

>AB=—^•(4a-3L) 
^ 12E1 


Solving for the stiffness: 


W _ 12 E1 lb 

>AB B^(3L-4a) 


For; 

UoDer Tier 


12 El , 

(Convert inches to feet) 

KB, ^-12 


a , (3 1, - 4 a j) 


L ower Tier 

12 El , 

KB2 =-=- 12 

82 (^'1^2” 


(per beam) 


KB 1 =2.947-10* — 
‘ ft 


KB 9 =2.947-10* 

ft 
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Input Number of Beams Per Tier; 


Beams | = 2 


Beams 2 = ^ 


K j = KB ] Beams j 


Tier Stiffness: K,-5 895-10* - 

* ft 


K 2 = KB 2 Beams 2 

K, = 1.179-10'^ - 

^ ft 


B. SYSTEM MASSES: (Includes Effective Sjirina Mass 
Upper Tier Wieght: W j = 2500 lbs Lower Tier Weight: W 2 = 2500 lbs 


* IT 


M I =77.64 lb 


M-, =77.64 lb 


Mass Ratio a: 


C. TIER NATURAL FREQUENCIES: 


.ower Tier 


f, 


1 K, 


2'rc JM 


D. SYSTEM NATURAL FREQUENCIES: 


f - J_. ^ 

^ 2n^M2 


fj =43.855 Hz f 2 =62.02 Hz 


fnl ^ = -F'J(‘'i)^^“ (^0^^(‘’2)^-4(fi)^+-“(fi;^+ (^ 2 )^] - 

V2 


f„j = 33.565 Hz 


fn2 ^■7’-](^l)^+“ (^ 2 )^+^* (•'l .'^-^ “ yf l)^+(^2)^1 - vh'j' 

'J2 


fjj 2 =81.033 Hz 



APPENDIX C3 


Mathcad Worksheet to Find Two DOF Coupled System Natural 
Frequencies as a Function of the Uncoupled Tier Frequencies 


Pick a Value For the Lower Tier Uncoupled Natural Frequency (F2); F2 - 48 Hz 
Specify the Mass Ratio; a = 1 0 


Solve for the Coupled System Natural Frequencies as a Function of F1: 

k ^ I . 100 Flj. =k 


FNl. 


r 


/t 


iiFl.i 

V , k- 


a 


■FI ‘ 
I k 


(F2r 


iFl 


s2 


a FI. 


(F2)- 


FI •F2'',^ 


FN2j^ =■ 


,/T 


JiFl,, 
'V ' k, 


\2 


a (Fl^;^+ (F2r-j[(Fg^+a ;F1^^-(F2)^ 


-i2Fl,F2V 


i =40.60 


PICK 

UPPER TIER COUPLED 


FI 

t 

12 

41 


12 

44 

45 

46 

47 

48 

12 

E 

51 

52 

53 

54 

22 

56 

57 

22 

59 


FNl FN2. 

_I. _i 


27.936 

68.728 

28 198 

69.792 

28.446 

70.872 

28.679 

71.969 

28.899 

73.081 

29 108 

74.207 

29.304 

75.348 

2949 

76.501 

29.666 

77.666 

29.832 

78.842 

29.989 

80.029 

30.138 

81.226 

30.279 

82.433 

30.413 

83.648 

30.54 

84.872 

30.661 

86.104 

30.775 

87.3'44 

30884 

88 59 

30.987 

89.843 

31.086 

91.103 

31.18 

92.368 


COUPLED SYSTEM NATURAL FREQUENCIES 

a = 1 F2 = 48 



UPPER TIER NATURAL FREQUENCY 
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APPENDIX C4 


MathCad Worksheet To Find Three DOF Fixture Ucoupled Frequencies: 


Mass Ratio: 


Input the Desired Coupled Natural Frequencies: (Hz) 

f„, ^30.00 -80 00 

Input The System Tier Weights: (lbs) 


1 =2500 lbs 

W 2 =2500 lbs 

W 3 =4500 

Wi 

Wo 

W , 

^ 1 

P ^ — 

Y = - 



W 3 

a = 1 

3 = 0 556 

Y = 0 556 


Supply the Program With An Initial Guess at f-i and ( 2 : 


f; ^40 Hz f 2 =40 Hz 


Given the initial guess for fi and ( 2 , the following equations are iterativley solved for the required 
values of f-| and f2 








l \ 


ff,)^-(l + a)+(f2)'(U(3) 




a)+ (f?; (1 + P) 




Y+ 1) 


\f2 


= Findff 1^2! 


Required Uncoupled Tier Natural Frequencies: f j = 54 077 Hz 


f 2 = 30.545 Hz 
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APPENDIX C5 


Mathcad Worksheet to Determine System Natural Frequencies Of Three DOF Fixture: 


A. SPRING STIFFNESS: 

The spring stiffness is based on the following beam model; 

y 




Simple Supports 

' Symetric Loads 


Supply the following Dimensions: 







1) Dimensions: 

Upi 

aer Tier 


Lower Tier 

Length: 

Ll 

= 50 

in 


L2 ^50 

in 

Load Application Point: 


= 13 

in 


82 ^13 

in 

Area Moment of Inertia: 

h 

= 11.3 

in 


I 2 =113 

4 

in 

Young's Modulus: 

E : 

= 30 10* 

psi 


E =30 10’ 

* psi 

2) Stiffness: 







From the beam bending 

equation: 

> AB 

6 El ' 

x^+ 3 a^- 3 l) 






_ W 

(4 a- 3 L) 


Let x=a and F=W/2; 



>■ AB 

12 E 1 





.. W 

12E 1 

lb 

Solving for the stiffness: 



K=- 

y 

AB 

(3 L- 4 a) 

in 


(Convert inches to feet) 


Upper Tier 


a 1 ^ (3 L 1 - 4 a , ) 


KB, ::- 


Lower Tier 

12EI, 


^ ® 2 ' 


(per beam) 


KB , =2 947-10*- 
‘ ft 


KB 2 = 2 947-10* 


I 






Upper Ti er 

Input Number of Beams Per Tier: Beams j =2 


Tier Stiffness: 


K ] =KB ) Beams j 
K , =5.895-iO'^ - 

' n 


B. SYSTEM MASSES: (Includes Effective Spring Mass) 
Upper Tier Weight: (W-j) Intermediate Weight: 


^'1 

= 2500 

W. =2500 



w. 

M , 

M. =—^ 

1 

32.2 

-- 32 2 

M , 

= 77.64 Ib-iL 

M . =77.64 


sec 


sec 


C. TIER NATURAL FREQUENCIES: Upper Tier 


IK 


f 


1 •; 


2 rt j 

f, =43.855 Hz 


D. System Natural Frequencies: 

Defining B as. 


Lower Tier 


Beams 2 -^ 


K 2 = KB 2 Beams 2 

K-, = 1179-10'^ - 

ft 

Anvil Table Weight: (W->) 

W 3 =4500 


M 


32.2 


M 3 = 139.752 lb 


Lower Tier 


ft 

scc‘ 




IK. 


2 jt ^jM. 


f 2 =62.02 Hz 


K , Kj K 1 +K. 
B =—+ 


M I M3 M2 


f. 


f 


:-L. l.|B- 


2 «^| 2 [ ^ 

First Frequency: fjjl=44 64 Hz 


2 4K,K2 ,1 

^ M M (Mi^M2^M3)I 

M ] M 2 M y ^ I 


*n 2 1:^ 

2n :. 


B + 


B- 


4.K,K2 , 


M ] M 2 M 3 


^M ] + M 2 -t- M 3 ': 


Second Frequency: 1^=88 528 Hz 
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APPENDIX C6 


Mathcad Worksheet to Find Three DOF Coupled System Natural 
Frequencies as a Function o{ the Uncoupled Tier Frequencies 

Pick a Value For the Lower Tier Uncoupled Natural Frequency (f 2 ) f o = 48 Hz 


Input The System Tier Wieghts 

W 1 = 2500 lb 

V '2 =2500 lb 

V,':, =4500 

and calculate mass ratios; 



V-’l 


a = — 

P =- 

Y "- 


W 2 


^'3 


a = 1 

P=0 556 

Y = 0 556 


Solve for the Coupled System Natural Frequencies as a Function of ri: 

Define Mass/Stiffness Constant: B. -,fi (I + a) * if'>’ 1 ^ (I + p) 


-k 






- 4 


I'f I 1 

\ k/ 


(pTy+ 1) 


I 


3 il |r , 4 ./f, ^-.(f2i’(P+Y* l)j 

-J2 [ ^ \ ‘k.' ^ ■ 


P/CK i - 30.. 50 

UPPER TiER COUPLED 


n. ^nl. ^ti 2 . 

I I I 


30 


31.593 

66.226 

31 


32.383 

66763 

32 


33 147 

67 33 

33 


33.882 

67.926 

34 


34.59 

68 553 

35 


35 269 

69.211 

36 


35.92 

69 898 

37 


36.542 

70 616 

38 


37 136 

71 365 

39 


21102 

72 142 

40 


38.242 

72949 



38 75J 

12 784 

42 


39.24 

74647 

43 


39.702 

75.537 

44 


40 138 

76452 

45 


40.552 

77.392 

46 


40.943 

78.356 

47 


■3’ 313 

79.342 

48 

y-o'.l 

30.35 r 

49 

■] ^^2 

SI.379 

50 

|4" .06 

32.427 


COUPLED SYSTE?JI NATURAL FREQUENCIES 

a = l f2=48 






150, 

y 



FN? 

■ 







— 








FN1 






_Q 


0 20 40 60 *0 100 


UPPER TIER NATURAL FREQUENCY 
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APPENDIX D. THREE DOF TUNED FIXTURE PROGRAM 


% Program THP.ELDOF.M 
% Lt. David M. Cox 
% Naval Postgraduate Sch 'Ol 


% This program plots the relative displacement,relative 
%velocity and absolute acceleration of a three DOF mass, 
%spring, damper system subjected to base excitation (half-sine 
%pulse). Additionally, the shock spectra of an undamped SDOF 
%system is calculated using the upper tier acceleration as a 
%base excitation to the SDOF system. 


%User Defined Variables 


HH=5 % 
Wtl=3000 % 
Wt2=2500 % 
Wt3=4500 % 
kbean=2.947e6 % 
kl=3*kbeam % 
k2=4*kbeam % 


2 eta=.04 


g=32.2 


% 


Drop height in feet 
VJeight of Upper Tier in Lbs 
Weight of Lower Tier in Lbs 
Weight of Anvil Table in Lbs 
Stiffness of one Beam Lbs/ft 
Upper Tier Stiffness in Lbs/ft 
Lower Tier Stiffness in Lbs/ft 
Critical Damping coefficient 
acceleration of gravity 


Base 


Acceleration 
% time 


%Fonnulate the 
dt=0.-0001; 

N=5000; 

Tblvel=sqrt(2*g*HH); 

PkAccel = ll.3 5*Tblvel 
T=N*dt; 

t=[0:dt:(N-1)*dt]; 
w=1000*pi; 

for i=l:N % loop to produce half-sine wave pulse 

if t(:.)< = 0.001 
Force(i) - (Wt3) 
else 

Force(i)= 0; 


% 

% 

% 

% 

% 


step for sampling 
Number of steps 
Velocity of Table after 
From Figure 19 
Total time record 
Time vector from t=0 to 


impact 


t=T 


'PkAccel*sin(w*t{i)) 


end 

end 


%Calculate the mass and mass ratio 
ml=Wtl/32.2 ; 
m2=Wt2/32.2 ; 
m3=Wt3/32.2 : 
m= [ ml, 0,0 ; 0 , in2,0 ; 0,0 , m3 ] ; 
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alpha=ml/m2; 
beta=m2 /it 3 ; 
gainma=ml/m3 ; 

%Calculate Circular Tier Natural Frequencies 
wl=sqrt(kl/rT\l); 
w2=sqrt (k2/in2) ; 

%Systein Natural Frequencies 

b=kl/iul + k2/m3 + (kl + k2 )/iti 2 ; 
wnl = 0 

wn2=sqrt (0.5* {b-sqrt (b''2- (4*kl*k2) * (ml+m2+m3 ) / (inl*m2*ni3 ) ) ) ) 
wn3=sqrt (0.5* (b+sqrt (b^2- (4*kl*k2) * (ml+in2+m3 ) / (inl*in2*m3 ) ) ) ) 

%Uncouple The Equations 
ull=l.0; 
u21=1.0; 
u3l=l.0; 
ul2=1.0; 

u22 = (kl-ml*vm2 ^2 ) / kl ; 

u32=(-kl + (kl + k2 - m2 * vm2^2)*u22)/k2; 
ul3=1.0; 

u23=(kl-ml*wn3 ^2)/kl; 

u33=(-kl + (kl + k2 - m2 * wn3'‘2) *u23 )/k2 ; 

U=[ull,ul2,ul3;u21,u22,u23;u31,u32,u33] 

%Decouple the System 
M=U'*m*U 

%Stiffness Matrix 

k=(kl,-kl,0;-kl,kl+k2,-k2;0,-k2,k2]; 

K-U'*k*U 

%Force Coefficients 
FC=U'*(0;0;1] 

%Solve the Uncoupled Equations of Motion 

% {nddotl) + [2*zeta*vml^2 ] * {ndotl) + [wnl''2 ] * {nl} =FCll*For ce/Ml 
% {nddot2 } + (2*zeta*wn2^2 ) * {ndot2 } + [wn2''2 ] * {n2 } =FC21*Force/M2 
% {nddot3) + [2*zeta*wn3^2] *{ndot3)+ [wn3''2] * {n3 ) =FC31*Force/M3 

%State Space Matrix for Eqiaation #1 
Al=[0,1;-(wnl^2),-2*zeta*wnl]; 

B1=[0;FC(1:1>/M(l:l)]; 

Cl=[1,0;0,1;-(wnl^2),-2*zeta*wnl); 

D1=[0;0;FC(1;1)/M(l;l)); 

[yl]=lsim(Al,Bl,Cl,Dl,Force,t); 

%State Space Matrix for Equation #2 
A2= [0,1; - (wn2''2 ) , - 2*zeta*wn2] ; 

B2=(0;FC(2:2)/M(5:5)]; 

C2=[l, 0;0, l;-(wn2''2; ,-2*zeta*wi-2 ] ; 

D2=[0;0;FC(2:2)/M(5:5)]; 

[y2]=lsim(A2,B2,C2,D2,Force, t); 
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%State Space Matrix for Equation #3 
A3= [0,1;- (vm3 ''2 ) , -2*zeta*wn3 ] ; 
B3=[0;FC(3:3) /M{9:9)]; 

C3=[1,0;0,1;-(wn3 ^ 2 ),-2*zeta*wn3]; 
D3=[0;0;FC(3:3)/M(9;9)]; 

[y3]=lsim(A3,B3,C3,D3,Force,t); 

%Couple the Equations 


XIDIS 

= 

ull*yl(:,1) 

+ ul2*y2(:,1) 

+ ul3*y3(;,1); 


XIVEL 


ull*yl(: ,2) 

+ ul2*y2(: ,2) 

+ Ul3*y3(:,2); 


XlACC 

= 

(ull*yl(: ,3) 

+ ul2*y2(: , 3 ) 

+ ul3*y3(:,3) 

) /g; 

X2DIS 

— 

u21*yl(: ,1) 

+ u22*y2(: ,1) 

+ u23*y3(: , 1) ; 


X2VEL 

= 

u21*yl(: ,2) 

+ u22*y2(:,2) 

+ u23*y3(:,2) ; 


X2ACC 

= 

(u21*yl(:,3) 

u22*y2 ( : , 3 ) 

+ u23 *y3(;,3) 

) /g; 

X3DIS 


u31*yl(:,1) 

+ u32*y2(;,1) 

+ u33 *y3( ; ,1) ; 


X3VEL 


u31*yl(; ,2) 

+ u32*y2(:,2) 

+ U33*y3( :,2) ; 


X3ACC 

=: 

(u31*yl(:,3) 

+ u32 *y2(:,3) 

+ u33*y3(;,3))/g; 


plot(t,XlACC),grid 

xlabeK'TIME (sec)'),ylabel(’ACCELERATION (g''s)') 

title { 'MODELED THREE-DOF FIXTURE UPPER TIER ACCELERA'J.-iON ’ ) 

gtext('5 Foot Hammer Drop (200 g''s)') 

gtext('Zeta = 0.04') 

gtext('Upper Tier Weight = 3500 lb') 

gtext('Lower Tier Weight = 2700 lb') 

gtext('Anvil Weight = 4500 lb') 

meta fig33.met 

pause 

%Upper tier velocity 
plot(t,XIVEL),grid 

xlabeK'Time (sec)'),ylabel('Ve.ocity - ft/sec') 
title('PREDICTED 3-DOF MODEL UPPER TIER VELOCITY') 
gtext('Wl=43.5 lbs; W2=57 lbs; W3=109 lbs') 
gtext('Drop Height=6 in; Zeta=0.02') 
meta threedof.met 
pause 

%Relative Displacements 

plot(t,(X3DIS-X1DIS)*12,t,(X3DIS-X2DIS)*12),grid 

title('PREDICTED RELATIVE DISPLACEMENT WITH RESPECT TO ANVIL 

TABLE•) 

xlabeK'Time (sec)'),ylabel('Displacement (in)') 

gtext('Wl=43.5 lbs; W2=57 lbs; W3=109 lbs') 

gtext('Drop Height=6 in; Zeta=0.02') 

gtext('Upper Tier'),gtext('Lower Tier') 

meta threedof.met 

pause 
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%Calculate the shock spectra 
NF=120; 

SF=1; 

DF=1; 

Preq=[SF:DF:NF]; 
vm=2*pi*Freq; 

F=[0;-1] ; 

H=0; 

for i=l:NF 

E=[0,l;-(wn(i)^2),0]; 

GG=[-(wn(i)"2) ,0] ; 

[spec]=lsim{E,F,GG,H,XlACC,t); 
maxspec(i)=max(abs(spec)); 

end 

plot(Freq,maxspec),grid,xlabel{'FREQUENCY (Hz)') 
ylabeK'ACCELERATION (g''s)*) 

title ( 'MODELED THREE DOF FIXTURE ACCELERATION SHOCK SPECTRA’ ) 
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